Initial data for Einstein's equations with superposed gravitational waves 
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A method is presented to construct initial data for Einstein's equations as a superposition of 
a gravitational wave perturbation on an arbitrary stationary background spacetime. The method 
combines the conformal thin sandwich formalism with linear gravitational waves, and allows detailed 
control over characteristics of the superposed gravitational wave like shape, location and propagation 
direction. It is furthermore fully covariant with respect to spatial coordinate changes and allows for 
very large amplitude of the gravitational wave. 
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I. INTRODUCTION 

Vacuum spacetimes which are perturbed away from 
stationary solutions of Einstein's equations are interest- 
ing in many different aspects. Using perturbations of 
flat space, one can examine (critical) collapse to black 
holes 0, 12 1 or investigate nonlinear interaction between 
gravitational waves. Perturbed black holes are expected 
to be produced by astrophysical events like binary black 
hole coalescence. Detailed understanding of the behavior 
of perturbed black holes, including the nonlinear regime, 
will be important for analyzing data from gravitational 
wave detectors like GEO, LIGO, TAMA and VIRGO. 
Numerical evolutions are the only known avenue to an- 
alyze Einstein's equation in general, three-dimensional, 
nonlinear situations. Moreover, non-stationary space- 
times without black holes or with just one perturbed 
black hole provide important testbeds and benchmarks 
for numerical evolution codes in a computational setting 
much simpler than a full binary black hole evolution. 

Such numerical evolutions require initial data rep- 
resenting perturbed spacetimes. Historically, Brill 
waves Pj are the most widely used approach to con- 
struct perturbations of Minkowski space with such a non- 
stationary component (e.g. They are based 
on certain simplifying assumptions, and allow for varia- 
tions of the gravitational wave perturbation through a 
freely specifiable function, commonly called q. Brill's 
idea has also been generalized to three dimensions and 
to black hole spacetimes (e.g., |il3|li|ll|ll|ll). All 
of these authors continue to encode the perturbation in 
a function q. It appears that, generally, this function is 
chosen rather ad hoc, its purpose mainly being to per- 
turb the spacetime in some way. While every (nonzero) 
choice for q leads to a perturbed initial data set, it is 
not clear what properties the perturbation has, nor how 
to control these properties. Given that q is often chosen 
to be bell-shaped (e.g., a Gaussian), it seems likely that 
the resulting perturbation is some vaguely localized lump 
of energy, rather than, say, a coherently traveling wave. 
Part of the motivation to use Brill waves was certainly 



that they lead to fairly simple equations which are easy 
to solve numerically. Since elliptic solvers have matured 
considerably over the last years (e.g. 0, 0, 0, Il7|). 
computational complexity is no longer a serious issue, 
and one is free to look for more general approaches, with 
easier to interpret properties of the resulting initial data 
sets. 

An obvious starting point are linearized gravitational 
waves [T§| . which then are incorporated into the solu- 
tion of the initial value problem 0, |23, Here, 
we continue along this line of thought, and propose a 
conceptually very simple method, which combines linear 
gravitational waves with the conformal thin sandwich for- 
malism p^ . Our basic idea is to build the linear grav- 
itational wave into the free data for the conformal thin 
sandwich equations. The method allows superposition 
of an arbitrary linear gravitational wave onto an arbi- 
trary background spacetime. The constructed data sets 
retain, at least qualitatively, the properties of the under- 
lying linear gravitational wave. Thus, properties of the 
perturbation to be inserted in the initial data set can be 
controlled easily by selecting the appropriate underlying 
linear gravitational wave solution. 

In Sec. ^ we present the method and discuss its rela- 
tionship to previous work 0, 0, |23, 0| . As an illustra- 
tion we superpose, in Sec. IIIII quadrupolar gravitational 
waves on Minkowski space and on a Schwarzschild black 
hole. We close with a discussion in Sec. IIVI 



II. METHOD 

Emplo ying the usual 3-1-1 decomposition of Einstein's 
eauations |23l l23| . the spacetime metric is written as 

('^^ds^ = -N^dt^ + g^j{dx' + P'dt){dx^ + l3'-'dt), (1) 

where gij represents the spatial metric on t=const. hy- 
persurfaces, and N and /?' denote the lapse function and 
shift vector, respectively. The extrinsic curvature, Kij is 

defined hy K — — -i ± -Ci^'^g, where ^^^g represents the 
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spacetime metric, n the future-pointing unit normal to 
the hypersurface, and _L the projection operator into the 
hypersurface. Einstein's equations then spht into evolu- 
tion equations, 



{dt-Cp)K,j = N {R,j-2K,kK'', 



(2) 



and constraint equations 

(K'^ -g'^K) = 0. 



K,jK'^ = 0, 



(3) 

(4) 
(5) 



Here, is the covariant derivative compatible with g^j, 
/^represents the Lie-derivative, and Rij denotes the Ricci 
tensor of gij. Furthermore, R and K denote the traces 
of the Ricci tensor and the extrinsic curvature, respec- 
tively, and we have assumed vacuum. Initial data for 
Einstein's equations consists of {gij,K^^). The difhculty 
in constructing such data lies in the requirement that 
the Hamiltonian and momentum constraints, Eqs. Q 
and jSJ, must be satisfied. 

One widely used formalism for constructing initial data 
is the conformal thin sandwich approach [22 . i25ii] . It is 
based on two neighboring hypersurfaces, their conformal 
three-geometries, and the instantaneous time-derivative 
of the conformal three-geometry. One introduces a con- 
formal spatial metric ijij, related to the physical spatial 
metric by 



9ij = 



4- 



(6) 



where ip is called the conformal factor. To construct ini- 
tial data, one chooses the conformal metric gij , its time- 
derivative 



(7) 



as well as the trace of the extrinsic curvature K and the 
conformal lapse N — iIj~^N. We note that Uij must be 
traceless, Uijg^-' = 0. Having made these choices, the 
Hamiltonian and momentum constraints take the form 
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-V'^V'if = 0. (9) 
o 



Here, and R are the covariant derivative compatible 
with gij and the trace of the Ricci tensor of gij , respec- 
tively, L denotes the longitudinal operator. 



and A^^ is defined as 



2N\ J 



(10) 



(11) 



Equations ^ and are elliptic equations for the con- 
formal factor ip and the shift /3*. After solving these 
equations for ip and /3% the physical initial data is given 
by Eq. © and by 



K,. 



+ -g'^K. 



(12) 



Instead of specifying N as part of the free data one 
can also set dtK. It is well known that this leads to an 
elliptic condition for the lapse- function (e.g. [25ll2^l2^ '): 

y\N^p')-{N^/)[^R+^rK' + li^-'A,,A^^ 

^-^^dtK-(3''dkK). (13) 

The second ingredient into the construction of per- 
turbed initial data is linear gravitational waves. In lin- 
earized gravity 18], the spacetime metric is written as 



Ah„ 



(14) 



where 



is the Minkowski-metric, A 1 a. constant, 



and hfj_„ = 0{1) the linear gravitational wave. (We sep- 

for later convenience.) 



arate the amplitude A from h 
In transverse-traceless gauge U, 



is purely spatial. 



h'i_io = 'iom = 0, transverse with respect to Minkowski 
space, V^hij = 0, and traceless, rf^hij — Q. To first or- 
der in the amplitude A, Einstein's equations reduce to 



nK, = 0, 



(15) 



where □ is the Minkowski space d'Alambertian. The 3-1-1 
decomposition of the metric l|14|l in transverse-traceless 
gauge is 



gij — fij A hij 

if = 0, 
A^= 1, 



(16) 
(17) 
(18) 



where /.y denotes the flat spatial metric. From the evo- 
lution equation for , Eq. Q , we find the extrinsic cur- 
vature 



V ■ 



(19) 



The spacetime metric (|14|l satisfies Einstein's equa- 
tions to first order in A. Consequently, {gij,K'^^) from 
Eqs. (|16|) and (|19|) will satisfy the constraints to linear 
order in A. Since we intend to increase A to order unity, 
this is not sufRcient, and we must solve the constraint 
equations. Because the spatial metric, Eq. H16() and its 
time-derivative, Ahij , are known, it seems appropriate 
that this information be incorporated into the constraint- 
solve. 

In light of the conformal thin sandwich formalism, it 
seems obvious to use Eq. (|16|l as conformal metric, and 
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to base the time-derivative of the conformal metric on 
hij'. 

g^j ^ l,j+Ahij, (20) 

Uij = Ah.j - ^g.jg'^'Ahki- (21) 

The second term in (|21|l ensures that Uij is tracefree with 
respect to gij. Because hij and hij are traceless, Eq. 
suggests the choice 

= 0. (22) 
The free data is completed by setting 

dtK = 0. (23) 

While the free data Eqs. (pn|l - lP5|l were motivated by 
a small perturbation, they can be used equally well for 
large amplitudes A (as long as solutions can be found). 
Hence, by increasing A, one can obtain nonlinearly per- 
turbed initial data sets. 

In writing down Eqs. H22|l and (|23|l . we have neglected 
terms of order 0{A^) on the right hand sides which arise 
because hij is traceless with respect to the flat metric 
fij, but not with respect to the perturbed metric fij + 
Ahij. Linearized gravity cannot determine such higher 
order terms. Since nonlinearities of Einstein's equations 
arise at the same order, and these nonlinearities are not 
accounted for in /i^j, we see no advantage to including 
O(A^) terms in Eqs. (j^U and if^ . We have also chosen 
to use Eq. (|23|l as free data and include Eq. (|13|) as a 
fifth elliptic equation. An alternative is to set TV = 1, 
and to solve only the four equations JHJl and 10. Both 
alternatives are identical to linear order in A. 

Equations H2()|l - H23|l . which result in a perturbation 
of Minkowski space, can be generalized to curved back- 
grounds easily by replacing the flat metric by a curved 
metric: Let q?„ and K'^ be the 3-metric and mean cur- 
vature of an asymptotically flat, spatial slice through a 
stationary spacetime (for example flat space or a Kerr 
black hole). Solve the conformal thin sandwich equations 
©, ®, and lO with the free data 

9^3 ^ 9"j+Ah,j, (24) 

u,j = Ah,j -^g.jg^^Ahku (25) 

K = K°, (26) 
dtK = 0. (27) 

We consider a few limiting cases 

• For A = the free data reduce to gij = gij , K = 
K^^Uij = dtK — 0. In this case, the underlying 
stationary spacetime is a solution of the conformal 
thin sandwich equations. 

• For A <C 1 and the wave hij located in the asymp- 
totically flat region of the hypersurface, linear the- 
ory is valid. The properties of the perturbation in 
the initial data set will be precisely those of the 
underlying linear wave hij . 



• For large A we will have a nonlinearly perturbed 
spacetime, our primary interest. Due to the non- 
linearity of Einstein's equations, the properties of 
such a strongly perturbed spacetime will differ from 
the linear wave. However, we expect that the qual- 
itative properties are unchanged. 

For constructing perturbed initial data on a curved 
background, one can, of course, also use a gravitational 
wave hij which represents a linear wave on the hack- 
ground gij, rather than on flat space. In that case, the 
A 1 limit approaches the underlying linear wave even 
if the underlying wave is located in the strong field re- 
gion. Since construction of linear waves on curved back- 
grounds is more complicated than on flat space, the deci- 
sion whether this is necessary for a particular application 
will depend on how closely the perturbation must match 
an exact linear wave in the limit A 1. Superposition 
of a flat space linear wave at intermediate separations 
from a black hole, say, lOM or 20M, should result in 
a gravitational wave which predominantly, although not 
exactly, retains the properties of the linear wave, which 
may be sufficient for many applications. 

Finally, we remark that our approach is related to and 
generalizes work by Abrahams & Eva ns IT^ , Shibata 
& Nakamura j20j and Bonazzola et al. |2j . Abrahams & 
Evans 0, Il9j assume axisymmety, and set a certain com- 
ponent of the extrinsic curvature (namely K^g in spherical 
coordinates) equal to the value appropriate for the lin- 
ear wave. Then they solve the momentum constraints for 
the remaining components of Kij . This procedure singles 
out a preferred coordinate system, while our method is 
covariant with respect to spatial coordinate transforma- 
tions. 

Shibata & Nakamura [23| use the extrinsic curvature 
decomposition 0, |2^ to construct initial data, rather 
than the conformal thin sandwich equations. By choos- 
ing a maximal slice, K = 0, Hamiltonian and momen- 
tum equation decouple, and the momentum constraint 
is solved by the analytical (transverse tracefree) extrin- 
sic curvature determined from the underlying linearized 
wave. This is a very elegant approach, since only the 
Hamiltonian constraint remains to be solved for the con- 
formal factor; however, decoupling of the Hamiltonian 
and momentum constraint happens only for slices with 
constant K . In order to construct perturbations of black 
hole spacetimes with non-constant K, the coupled ellip- 
tic equations within the extrinsic curvature decomposi- 
tion have to be solved. This extension is conceptually 
straightforward and would lead to a method parallel to 
the one presented here, but using the extrinsic curvature 
decomposition rather than the conformal thin sandwich 
formalism. 

Bonazzola et al. 21j, finally, uses the conformal thin 
sandwich formalism to superpose a time-symmetric lin- 
earized wave {dtgij = 0, i.e. a superposition of an in- 
coming and an outgoing wave) on flat space. This work 
employs Dirac-gauge and relies heavily on spherical coor- 
dinate systems in the reduction of the remaining degrees 



4 



of freedom to two scalar functions. Our method, in con- 
trast, can be used with any spatial coordinates, which 
is of particular importance for superposition of gravita- 
tional waves on black hole backgrounds, which may not 
be available in Dirac-gauge. Furthermore, while we use 
Teukolsky waves below as an example, our method can 
also applied to different linearized waves, for example 
a spherical wave superposed with a plane wave, which 
would be more difficult to implement in the approach 
taken in Ref. 111. 



III. NUMERICAL RESULTS 

A. Quadrupole waves 

We illustrate the general method introduced in Sec. ^ 
with linearized quadrupole waves as given by Teukol- 
sky p^. This reference explicitly presents even parity 
waves, which are superpositions of / = 0, 2, and 4 modes, 
as well as odd parity waves, which are constructed as 
superpositions of ^ = 1 and 3. For each parity, there 
are five independent modes, corresponding to azimuthal 
quantum number m = ±2, ±1, 0. The even parity outgo- 
ing wave has a spacetime line-element 



(1 + Afrr)dr^ + {2Bfr0)rdrde 
{2Bfre)r sin edrd(p 

[2{A-2C)fg^] sin 0ded(j) 



(28) 



with radial dependence given by 
3^(1) 3Fl 



2 • 2 

r sm 



A = 3 
B = - 



f (3) ^ 3f (2) ^ Qp{l) ^ Qp 

F(4) 2i^(3) 9^(2) 21i^(i) 21i^' 



where 



d"F(a;) 



(29) 
(30) 
(31) 

(32) 



F{x) = F{t — r) describes the shape of the wave. The 

(2.) 

functions fm ■ ■ ■ , f^^ depend only on angles {9, (/>); they 
are given explicitly in Ref. for each azimuthal quan- 
tum number M . Ingoing quadrupole waves are obtained 
by replacing F(t — r) with a function oft + r, and revers- 
ing the signs in front of odd derivatives of F in Eq. (|32|l . 
Reference '2^ gives also the metric for odd parity waves. 
From Eq. H28|l . one can easily extract hij and hij. 




FIG. 1: Domain decomposition in R^. A cube covers the 
central region which is not covered by the spherical shells. 



B. Flat space with ingoing pulse 

We consider a perturbation of flat space, gfj — fij, 
— 0. We choose the even parity, to = ingoing 
mode. The shape of the pulse is taken as a Gaussian 



F{x) 



(33) 



of width w = 1 and with an initial radius of xq = 20. 

Equations JHl, Q, and ifT^ are solved with the pseu- 
dospectral elliptic solver described in |l(i| |. The domain 
decomposition used in the elliptic solver is shown in Fig- 
ure ^ We use three spherical shells with boundaries at 
radii r = 1.5,16,24, and 10^, so that the middle shell 
is centered on the gravitational wave. The inner shell 
does not extend to the origin, since the regularity condi- 
tions at the origin of a sphere are not implemented in the 
code. Instead, we place a cube on the origin which over- 
laps the innermost spherical shell. The solutions of the 
constraint equations turn out to be very smooth except 
for high frequency radial variations at the location of the 
gravitational wave pulse, cf. Eq. (|33|l . Therefore, the 
accuracy is completely determined by the number TV^ of 
radial basis functions in the middle spherical shell, which 
is chosen significantly larger than the number of basis 
functions in the other subdomains. At the highest res- 
olution, there are — 84 radial basis functions in the 
middle shell, but only 28 in the other two shells. Fur- 
thermore, the angular resolution of all shells is i = 15 
and the cube has 18 basis functions in each dimension. 

Figure|21presents the residuals of Hamiltonian and mo- 
mentum constraints, Eqs. Q and |SJ for the linear grav- 
itational wave without solving the constraints, i.e. upon 
direct substitution of Eqs. (|16|1 and H19I) into the con- 
straint equations. As expected, the residual is 0{A^), 
confirming that the quadrupole wave is indeed a solution 
of linearized gravity. 

We now solve the conformal thin sandwich equations 
with the free data (|24|I - H27() for different A, and compute 
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FIG. 2: Constraint violation of linear gravitational wave in 
flat background prior to solving the constraints. 
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FIG. 3: Convergence of the elliptic solver for different ampli- 
tudes A. Plotted is the residual in the Hamiltonian constraint 
(root mean square) vs. the number of radial basis functions 
in the middle spherical shell. 



the ADM energy for each solution, 



^^ADM — [ {gij.j - 9j3.i) if' Si 



(34) 



Figures 13 and 01 plot the residual of the Hamiltonian con- 
straint and the error of the ADM energy versus iV^. Ex- 
ponential convergence is apparent, until roundoff limit is 
encountered around ~ 10"-'^''. The seemingly large value 
of this number (when compared to the usual double preci- 
sion floating point accuracy of ~ 10~^^) is a consequence 
of the many numerically computed derivatives that en- 



FIG. 4: Convergence of the elliptic solver for different ampli- 
tudes A. Plotted is the difference of the ADM energy to the 
next higher resolution solution vs. the number of radial basis 
functions in the middle spherical shell. 



1 - 




0.3 



FIG. 5: ADM energy of an ingoing Gaussian pulse in flat 
space. The dashed line indicates the low-amplitude quadratic 
behavior. 



ter the calculation: To compute R in Eq. Q), second 
numerical derivatives of the conformal metric Eq. Hlt)|) 
are taken. After solution of the elliptic equations, the 
physical metric gij is assembled, and the Hamiltonian 
constraint Eq. Q is evaluated using second numerical 
derivatives of gij. Each of these numerical differentia- 
tions increases the roundoff error by a factor of order the 
number of basis functions. The increase in roundoff error 
with the number of basis functions can be clearly seen in 
Fig. 3 
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the origin with radius r w. Black holes usually form 
for systems with mass to size ratio of order unity; here, 
Eadm/w ~ 2.8, so that black hole formation appears 
very likely once the pulse is concentrated at the origin. 

These data sets could be used to examine critical col- 
lapse to a black hole, repeating Abrahams & Evans 
and extending it to genuinely three-dimensional collapse 
by choosing m 7^ in the underlying quadrupole wave. 
These datasets also provide a testbed for evolution codes 
in situations where the topology of the horizons changes. 



C. Black hole with gravitational wave 

As a second example of the flexibility of our method, 
we superpose a gravitational wave on a black hole back- 
ground. The background spatial metric and trace of the 
extrinsic curvature are set to a Schwarzschild black hole 
in Eddington-Finkelstein coordinates. 



FIG. 6: Cuts through the equatorial plane of the A — 0.3 data 
set of Fig. El The large plot shows lapse and conformal factor, 
the inset shows the scalar curvature of the 3-metric. 
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-3/2 



(35) 
(36) 



We now turn our attention from the convergence prop- 
erties to the actual solutions of the constraint equations. 
For small A, we find that -0 — 1 is proportional to 
everywhere. This is expected, because ip — I corrects the 
conformal metric to satisfy the Hamiltonian constraint. 
As the constraint violation is proportional to A^, so is 
this correction. Fig. [S] proves that one can clearly achieve 
initial data sets with a significant energy content. At 
low amplitudes, i^ADM is proportional to A^, as one ex- 
pects given that ip — 1 is proportional to A^. At high 
amplitudes, however, ii^ADM grows faster than A^, indi- 
cating that the non-linear regime with self-interaction is 
reached. For A > 0.3, the elliptic solver fails to converge. 

We now discuss the data set with amplitude, A — 0.3 
in more detail. Its ADM energy is i?ADM = 2.858. Fig- 
ure El presents cuts through the conformal factor ip, lapse 
iV and the scalar curvature of the physical 3-metric, 
^■^••i?. Conformal factor and lapse deviate significantly 
from unity confirming that the solution is indeed deep 
in the nonlinear regime. The scalar curvature is virtu- 
ally zero everywhere except within a spherical shell with 
18 ^ r < 20. Although the linearized wave is based on 
a Gaussian profile Eq. I|33|l . even finer features are in- 
troduced into the linearized wave (and into the scalar 
curvature plotted in Fig.O because of the derivatives in 
Eqs. (|29|I - H31() . Resolution of these fine features necessi- 
tates the high radial resolution of Nr. 

The gravitational wave is concentrated in a spherical 
shell of width w = 1. The underlying linear wave is 
purely ingoing, so it seems reasonable that the gravi- 
tational perturbation in the physical, nonlinear space- 
time is also predominantly ingoing. Neglecting disper- 
sion, the wave will concentrate in a sphere centered at 



where — x'' /r, and ~ SijX^x^ . 

We choose an odd, ingoing m — quadrupole wave 
with Gaussian shape, Eq. (|33|l at location Xq = 15 and 
width w = 1. The metric is singular at the origin, there- 
fore we excise at an inner radius of 1.5M (which is inside 
the horizon). At this inner boundary, we impose simple 
Dirichlet boundary conditions appropriate for the unper- 
turbed black hole: ip = I, and N — Nq and (3^ = (3q^ with 
lapse and shift for Eddington-Finkelstein given by 



-1/2 



r / 



(37) 
(38) 



Perturbed initial data sets are constructed for vari- 
ous values of A, and Fig. [71 demonstrates convergence 
of the solutions. In this case, the resolution is deter- 
mined by two factors, namely how well the gravitational 
wave is resolved in the middle shell, and how well the in- 
ner shell resolves the background solution (g? , K^^ ) . For 
perturbations of Minkowski space, the latter was trivial 
(any expansion resolves the constant Minkowski back- 
ground), however, here it is the limiting factor for small 
A or low resolutions, whereas for large A and high reso- 
lutions, both effects are about equally important. At the 
highest resolution, the number of radial basis functions 
in each shell is (from inner to outer) 48, 84 and 28, and 
the angular resolution is unchanged from before, L = 15. 
In each resulting initial data set, the apparent horizon 
is located with the a ppa rent horizon finder implemented 
and tested in |2^ IsollSlj , and the apparent horizon mass 
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FIG. 7: Black hole with superposed gravitational wave: 
Residual of the Hamiltonian constraint vs. radial number 
of basis functions in middle spherical shell. 
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FIG. 8; Black hole with superposed gravitational wave. 



is computed from the area of the apparent horizon via 



A 



AH 



(39) 



Figure |H1 presents the ADM energy and the apparent 
horizon mass of the central black hole as a function of 
the amplitude of the gravitational wave. The apparent 
horizon mass is fairly independent of A indicating that 
the horizon of the central black hole is only slightly per- 
turbed by the gravitational wave. However, the ADM 
energy, which measures the total energy in the hypersur- 
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FIG. 9: Apparent horizon mass during an evolution of a per- 
turbed black hole spacetime. The dashed line indicates i?ADM 
as computed from the initial data set. 



face, depends strongly on A; for large A, 



Mah 



>2.5, 



(40) 



indicating that a significant amount of gravitational en- 
ergy resides outside the black hole. 

To support our assertion that the superposed initial 
data set retains the features of the underlying gravita- 
tional wave, we present a preliminary evolution of a black 
hole with superposed ingoing gravitational wave [3^ . 
The initial data for the evolution is identical to the 
data sets used in Figure |H1 with the one exception that 
the gravitational wave has even parity. Figure El shows 
the apparent horizon mass as a function of evolution 
time. All quantities are scaled such that the unperturbed 
Schwarzschild black hole has unit mass. For t < 10, 
A/ah is constant, its value being that from the initial 
data set. Between 12 < t < 15, A/ah increases rapidly 
to an asymptotic value of A/^jj « 1.084. The ADM en- 
ergy of the initial data set was /?adm = 1.0845. Appar- 
ently, the ingoing gravitational wave outside the black 
hole falls into it, increasing the area of the apparent hori- 
zon. The final apparent horizon mass is very close to the 
ADM energy, and the growth of Mah happens during a 
time-interval comparable to the width of the initial pulse. 
Thus it appears that a large fraction of the wave is co- 
herently ingoing and falls into the black hole. 
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IV. DISCUSSION 

We propose a conceptually very clear method to 
construct spacetimes containing gravitational radiation 
which combines the conformal thin sandwich formalism 
with linear gravitational waves. For small amplitudes, 
the gravitational perturbation in the resulting initial data 
sets retains the characteristic features of the underlying 
linear wave, allowing for easy control of the properties of 
the gravitational wave perturbation. For strong ampli- 
tudes, nonlinearities of Einstein's equations are impor- 
tant, but we expect that the solutions still retain quali- 
tatively the properties of the underlying linear wave. 

To illustrate the method, we superpose quadrupolar 
gravitational waves onto Minkowski space, and onto a 
Schwarzschild black hole. In both cases, initial data with 
a large amount of gravitational energy in the perturba- 
tion can be constructed. 

Numerically, these initial data sets provide test-beds 
of evolution codes in situations away from stationarity. 
The mass of a central black hole changes -it may even 
double- when a large gravitational wave falls into it; can 
current gauge conditions handle this situation? If a grav- 
itational wave collapses to a black hole, horizons appear, 
and evolution codes using black hole excision must ac- 
commodate this change. Furthermore, spacetimes with 
outgoing gravitational wave perturbations are ideal test- 
beds for gravitational wave extraction algorithms, or con- 



straint preserving boundary conditions [s^ . 

Physically, ingoing gravitational wave pulses in 
Minkowski space, like the ones presented in Figs. |S1 and El 
could be used to examine critical collapse, including 
the genuinely three-dimensional regime with m 7^ 0. 
The black hole initial data sets with ingoing gravita- 
tional wave pulses (cf. Fig. |SJ) would be useful to ex- 
amine scatterin g of the gravitational wave at the black 
hole |3J, [Sa, [Sy, |33 : What fraction of the gravitational 
wave is scattered and reaches infinity? Which multipole 
moments are excited in this process? This example can 
also be generalized to spinning black holes, off-centered 
gravitational waves, or gravitational waves with m 7^ 0. 
Interesting questions in these scenarios would include, 
whether one can impart linear or angular momentum on 
the black hole. 
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